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Abstract. Fractional calculus and q-deformed Lie algebras are closely related. 
Both concepts expand the scope of standard Lie algebras to describe generalized 
symmetries. For the fractional harmonic oscillator, the corresponding q-numbcr 
is derived. It is shown, that the resulting energy spectrum is an appropriate tool 
e.g. to describe the ground state spectra of even-even nuclei. In addition, the 
equivalence of rotational and vibrational spectra for fractional q-deformed Lie 
algebras is shown and the B a (E2) values for the fractional q-deformed symmetric 
rotor are calculated. 
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1. Introduction 

The combination of concepts and methods developed in different branches of physics, 
has always led to new insights and improvements. In nuclear physics for example, 
the description of rotational and vibrational nuclear spectra has undoubtedly been 
influenced by concepts, which were first developed for molecular spectra. 

The intension of this paper is similar. We will show, that the concept of q- 
deformed Lie algebras and the methods developed in fractional calculus are strongly 
related and may be combined leading to a new class of fractional deformed Lie algebras. 

For that purpose, we will first introduce the necessary tools for a fractional 
extension of the usual q-deformed Lie algebras and present the fractional analogue 
of the standard q-deformed oscillator. 

As an application we will use the fractional q-deformed oscillator for a description 
of the low energy excitation ground state band spectra of even-even nuclei. 

In addition, the equivalence of rotational and vibrational spectra for the fractional 
q-deformed Lie algebras is demonstrated and the B a (E2) values for the fractional q- 
deformed symmetric rotor are discussed. 

2. q-deformed Lie algebras 

q-deformed Lie algebras are extended versions of the usual Lie algebras pQ. They 
provide us with an extended set of symmetries and therefore allow the description of 
physical phenomena, which are beyond the scope of usual Lie algebras. In order to 
descibe a q-deformed Lie algebra, we introduce a parameter q and define a mapping 
of ordinary numbers x to q-numbers e.g. via: 




q x - q 



— X 



(1) 



q - 



-i 
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which in the limit q — > 1 yields the ordinary numbers x 

lim [x] = x (2) 

q^l » 

Based on q-numbers a q-derivative may be defined via: 

With this definition for a function /(x) = x™ we get 

D« = a;"- 1 (4) 

As an example for q-deformed Lie algebras we introduce the q-deformed harmonic 
oscillator. The creation and anihilation operators a', a and the number operator N 
generate the algebra: 

[N, a 1 "] = a f (5) 

[N,a] = - a (6) 

aa) - ^Va = q ±N (7) 

with the definition ([T]) equation (J7J) may be rewitten as 

t„ - 



a 1 a 



t _ 



M, (8) 
[N + l] q (9) 

Defining a vacuum state with a|0 >= 0, the action of the operators {a, a', N} on the 
basis \n > of a Fock space, defined by a repeated action of the creation operator on 
the vacuum state, is given by: 

N\n > =n\n> (10) 
a f |n > = 0n+ l] q \n + 1 > (11) 

a|" > = \/N^I" - 1 > (12) 
The Hamiltonian of the q-deformed harmonic oscillator is defined as 

H = —(aat+a^a) (13) 

and the eigenvalues on the basis |n > result as 

E q (n) = ^-([n] q + [n + l} q ) (14) 

3. Fractional calculus 

The fractional calculus [2], [5] provides a set of axioms and methods to extend the 
coordinate and corresponding derivative definitions in a reasonable way from integer 
order to arbitrary order a, where a is a real number: 

<*5>~ <*■■£> <15) 

The definition of the fractional order derivative is not unique, several definitions e.g. 
the Riemann, Caputo, Weyl, Riesz, Feller, Griinwald fractional derivative definition 
coexist 
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We will apply the Caputo derivative 

1 /" .. U 



rde(x-0~ a 4/(0 forO<a<l 



r(i - a) J 

1 r, 

r(2 - a) A 
For a function set f(x) — x na we obtain: 

{rfl + na) („ ,i„ 
— v . {--x {n - 1)a for n > 
r(l + (n - I)*) (17) 
for n = 

Let us now interpret the fractional derivative parameter a as a deformation parameter 
in the sense of q-deformed Lie algebras. Setting \n >= x na we define: 

( r(l+na) . 

T ,., fr- n> for n > 

H>>= < r(l + (n-l)a)' (18) 

[ for n = 

Indeed it follows 

lim [n] — n (19) 

a — >1 

The more or less abstract q-number is now interpreted within the mathematical 
context of fractional calculus as the fractional derivative parameter a with a well 
comprehended meaning. 

The definition (|18[) looks just like one more definition for a q- deformation. But 
there is a significant difference which makes the definition based on fractional calculus 
unique. 

Standard q-numbers are defined more or less hcuristically. There exists no 
physical or mathematical framework, which determines their explicit structure. 
Consequently, many different definitions have been proposed in the literature (see 
e.g. P). 

In contrast to this diversity, the q-deformation based on the definition of the 
fractional derivative is uniquely determined, once a set of basis functions is given. 

As an example we will derive the q-numbers for the fractional harmonic oscillator 
in the next section. 



4. The fractional q-deformed harmonic oscillator 

The transition from classical mechanics to quantum mechanics may be performed by 
canonical quantization [TD],[TT]. The classical set of canonical observables {x,p} is 
replaced by a set of quantum mechanical observables {x,p}. In fractional calculus, 
the space coordinate representation of these operators is given by [12] : 

(l-a) 

x a (20) 



p= -i I — ) mcD% (21) 



mc 



The attached factors (h/mc)^ 1 and (h/mc) a mc respectively ensure correct length 
and momentum units. 
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With these operators, the classical Hamilton function for the harmonic oscillator 

ffciass = + -muj 2 x 2 (22) 
2m 2 

is quantized. This yields the Hamiltonian H a 

Ha = & i + \™*# W 
and the corresponding stationary Schrodinger equation is explicitely given by: 

- (-s (£)° mVlKD ° + r" 2 - £ * (24) 

Introducing the variable £ and the scaled energy E 1 : 

£ = 7ia;.E' (26) 

we obtain the stationary Schrodinger equation for the fractional harmonic oscillator 
in the canonical form 

H a ^n(0 = l(-Dl a + e a )^n(0 = £'(n,a)*„(0 (27) 

Laskin [13] has derived an approximate analytic solution within the framework of the 
WKB-approximation, which has the advantage, that it is independent of the choice of 
a specific definition of the fractional derivative. We adopt his result: 

E\n,a) = {\+n) a ^(^^ a n = 0,1,2,... (28) 

In view of q-deformed Lie-algebras, we can use this analytic result to derive the 
corresponding q-number. With l|14p the q-number is determined by the recursion 
relation: 

[n + l} a = 2E'(n,a)-[n] a (29) 

The obvious choice [0] =0 for the initial condition leads to an oscillatory behaviour 
for [n] for a < 1. If we require a monotonically increasing behaviour of [n] for 
increasing n, an adequate choice for the initial condition is 

[0]„ = 2 1+Q W 2 [-^) (C(-«, \) C(-a, \)) (30) 
The explicit solution is then given by: 

Ma = 2 1+a - Q/2 (^r^)" («-«. 3 + \) - C(-«, 5 + 1 )) (3D 

where £(s, ie) is the incomplete Riemann or Hurwitz zeta function, which is defined 
as: 

oo 

a-s,x) ( 32 ) 

fc=0 

and for a = m S 0, 1, 2, ... it is related to the Bernoulli polynomials B m via 

C(-m,a;) = --; — ^— —B m+1 (x) (33) 
(to + 1) 




Figure 1. Plot of the ground state wave function <I/q+(£) of the fractional 
harmonic oscillator j2Z}i solved with the Caputo fractional derivative l|16[ l for 
different a. 



Of course, the vacuum state |0 > is no more characterized by a vanishing expectation 
value of the annihilation operator, but is defined by a zero expectation value of the 
number operator, which is the inverse function of the fractional q- number (|31[) : 

iV|0>= ([n] a ) _1 |0>=n|0>=0 (34) 

Setting a = 1 leads to [ti] j = n and consequently, the standard harmonic oscillator 
energy spectrum E'(n,a = 1) = (1/2 + n) results. For a = 2 it follows 

£?'(„, a = 2)= -for (£jfj) {\+nf (35) 

87T 3 ,1 9, ,„„. 

= r7!^ ( 4 +n + n) (36) 

2tt 3 8tt 3 , , , N 

which matches, besides a non vanishing zero point energy contribution, a spectrum 
of rotational type E rot = 1(1 + 1), I — 0,1,2,... Unlike applications of ordinary 
q-numbers, this result is not restricted to a finite number n, but is valid for all n. 

In figure ([T]) we have plotted the ground state wave function of the fractional 
harmonic oscillator (|27[) . obtained numerically for the Caputo fractional derivative 
(|16[) for different a. While for a w 1 the wave function behaves like exp(— ^ 2 /2), in 
the vicinity of a w 2 it looks more like a Bessel type function. The fractional derivative 
coefficient a acts like an order parameter and allows for a smooth transition between 
the ideal rotational (a = 2) and ideal vibrational (a = 1) limits. 

The properties of the fractional harmonic oscillator are therefore well suited to 
describe e.g. the ground state band spectra of even-even nuclei. We define 

Ey h (a, m £ h ,a£ h )=m v j h + a v j h ([J] a + [J+l] a )/2 J = 0,2,4,... (38) 

where mj lb mainly acts as a counter term for the zero point energy and <2g lb is a 
measure for the level spacing and use (|38]> for a fit of the experimental ground state 
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Figure 2. The upper row shows energy levels of the ground state bands for 
156 Gd, 19B Pt, 110 Cd and 218 Po, which are plotted for increasing J. Squares 
indicate the experimental values given in |14 | -|17 | . The full line indicates the 
optimum fit for the fractional q-deformed harmonic oscillator using J38D and 
dashed lines give the best fit for the fractional q-deformed SU a (2) symmetric 
rotor according to 1421 . In the lower row the corresponding backbending plots 
are shown, which allow for the determination of the maximum angular momentum 
•/max for a valid fit below the onset of alignment effects. 



band spectra of 156 Gd, 196 Pt, 110 Cd and 218 Po which represent typical rotational-, 
7-unstable-, vibrational and linear type spectra [T4"]-[T7]. 

In table [1] the optimum parameter sets for a fit with the experimental data are 
listed. In figure [2] the results are plotted. Below the onset of microscopic alignment 
effects all spectra are described with an accuracy better than 2%. Therefore the 
fractional q-deformed harmonic oscillator indeed describes the full variety of ground 
state bands of even-even nuclei with remarkable accuracy. 

5. The fractional q-deformed symmetric rotor 

In the previous section we have derived the q-number associated with the fractional 
harmonic oscillator. Interpreting equation (131 [) as a formal definition, the Casimir- 
operator C(SU Q (2)) 

C(SU a (2)) = [J] a [J + l] a J = 0,1, 2,... (39) 

of the group SU Q (2) is determined. This group is generated by the operators J+, Jo 
and J_, satisfying the commutation relations: 

[Jo, J±] = ± J± (40) 
[J+,J_] = [2J ] Q (41) 
Consequently we are able to define the fractional q-deformed symmetric rotor as 

Er(a,ml \ a ^)= m r + ar[J} a [J+l} a J = 0,2,4,... (42) 

where m™' mainly acts as a counter term for the zero point energy and a ot is a 
measure for the level spacing. 

For a — 1, E lot reduces to E TOt = m ot + a ot J(J +1), which is the spectrum of 
a symmetric rigid rotor. For a = 1/2 we obtain: 

lim (E}°1 2 - E l j ot ) = a u ot 7 r/2 = const (43) 



q-deformed Lie algebras and fractional calculus 



7 



miclid 



92 Gd 6 4 

196p f 
118 Ft 78 

llOrVI 
134 P °84 



Table 1. Listed are the optimum parameter sets (a, fiQ lb , according 
for the fractional harmonic oscillator for different nuclids. The maximum angular 
momentum Jmax for a valid fit below the onset of alignment effects is given as 
well as the root mean square error AE between experimental ( |14| -[T7 ) ) and fitted 
energies in %. 



a 


a# h [keV] 


ml lh {keV] 


J max 


AE[%] 


1.795 


15.736 


-14.136 


14 


1.48 


1.436 


91.556 


-39.832 


10 


0.16 


1.331 


197.119 


-87.416 


6 





0.801 


357.493 


-193.868 


8 


0.06 



which is the spectrum of a harmonic oscillator. We define the ratios 

( cwib pvib \ 

= TJL gU (44) 



«vib 


/ JTwib 

V h j 


- £ vib ) 


J, a 




- ^o vib ) 


■rot 


( TTTOt 


- ^5 ot ) 


J, a — 


( ZTTOt 


- £o ot ) 



E>rot — v •> U I (AK\ 

±L J,a (_B| ot — E rot ) 

which only depend on J and a. A Taylor series expansion at J = 2 and a — 1 leads 
to: 



K% = 1 + {J - 2)(0.10 - 0.04(J - 2)) + {a - 1)(J - 2)(0.101 - 0.020(J - 2)) (46) 

RT j%/2 = 1 + (J-2)(0.11 - 0.04(J-2)) + {a - 1)( J - 2)(0.095 - 0.016( J - 2)) (47) 

A comparison of these series leads to the remarkable result 

R$ +0(J\a 2 ) (48) 

Therefore the fractional q-deformed harmonic oscillator (f3"5)) and the fractional q- 
deformed symmetric rotor (|42p generate similar spectra. As a consequence, a fit of 
the experimental ground state band spectra of even-even nuclei with (|3S[) and (|4"2"]) 
respectively leads to similar results, as demonstrated in figure[2j There is no difference 
between rotations and vibrations any more, the corresponding spectra are mutually 
connected via relation (|4*5j) . 

Finally we may consider the behaviour of B (E2)-v&\\ies for the fractional q- 
deformed symmetric rotor. Using the formal equivalence with q-deformation, these 
values are given by [18] : 

B a (E2, J + 2 -> J) - ^o [2]j2J + 2]j2J + 3]a[2J + 4]Q[2J + 5]a (49) 

In figure[3]these values are plotted, normalized with 100/BQ,(i?2; 2 + — » + ). Obviously 
there is an saturation effect for increasing J. 

6. Conclusion 

We have shown, that fractional calculus and q-deformed Lie algebras are closely 
related. Both concepts expand the scope of standard Lie algebras to describe 
generalized symmetries. While a standard q-deformation is merely introduced 
heuristically and in general no underlying mathematical motivation is given, the q- 
deformation based on the definition of the fractional derivative is uniquely determined, 
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Figure 3. B(E2)-values for the fractional q-deformed SU a (2) symmetric rotor 
according to J29), normalized with 100/B a (E2; 2+ -» 0+). 



once a set of basis functions is given. For the fractional harmonic oscillator, we derived 
the corresponding q-number and demonstrated, that the resulting energy spectrum is 
an appropriate tool e.g. to describe the ground state spectra of even-even nuclei. 

In addition, the equivalence of rotational and vibrational spectra for fractional 
q-deformed Lie algebras has been demonstrated and we have calculated the B a (E2) 
values for the fractional q-deformed symmetric rotor. 

The results derived encourage further studies in this field. The application 
of fractional calculus to phenomena, which until now have been described using 
q-deformed Lie algebras will lead to a broader understanding of the underlying 
generalized symmetries. 
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